Measurements of Semileptonic Kl Decay Form Factors 
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We present new measurements of Kl semileptonic form factors using data collected in 1997 by the 
KTeV (E832) experiment at Fermilab. The measurements are based on 1.9 million Kl ~* w e^v 
and 1.5 million Kl — » 7r l^v decays. For f+(t), we measure both a linear and quadratic term: 
\' + = (20.64±1.75)xl0~ 3 and A'| = (3.20±0.69)xl(T 3 . For f (t), we find A = (13.72±1.31)xKT 3 . 
These form factors are consistent with K form factors, suggesting that isospin symmetry breaking 
effects are small. We use our measured values of the form factors to evaluate the decay phase space 
integrals, I e K = 0.15350 ± 0.00105 and 7j£ = 0.10165 ± 0.00080, where errors include uncertainties 
arising from the form factor parametrizations. 
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I. INTRODUCTION 



Measurements of semileptonic form factors provide 
unique information about the dynamics of strong inter- 



with precision comparable to the older experiments. The 
measurement presented in this paper provides another 
check of the older results with better precision. 



actions, 
and Kl 



The form factors for the decays Kl 



^e -1 



z l^v also are needed to determine the decay 



Jjr and 7^-, and thus are important 



phase space integrals 

for measuring the Cabibbo-Kobayashi-Maskawa [j], Q 
matrix element \V US \ 0- lj|- In this paper, we report a 
new measurement of the semileptonic form factors based 
on data collected by the KTeV experiment at Fermilab. 

The dependence of the semileptonic form factors on t, 
the four-momentum transfer to the leptons squared, was 
under serious experimental and theoretical study from 
the late 1960s to the early 1980s; the result of this effort 
was the precise measurement of the form factor in the 
decay Kl — > 7r ± e =F z/ ; 5j. Nearly twenty years later, this 
form factor measurement was confirmed by CPLEAR 6] 
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The form factors in the decay Kl —* tt ^v, on the 
other hand, have been a long standing subject of experi- 
mental controversy [jj . The inconsistency of different ex- 
perimental results leads to a large uncertainty in 1%, re- 
ducing the usefulness of the Kl — > 7r ± /i =F z/ partial width 
in the extraction of \V US \. The new KTeV measurement 
reduces the uncertainty in 7^- to the same level as in I e K . 
The improved precision of 1% allows a sensitive test of 
the consistency between the phase space integrals and 
the ratio of semileptonic partial widths TKfi,3/^Ke3- 

This paper is organized as follows. Section [B] intro- 
duces the phenomenology of semileptonic form factors. 
Section liTTl describes the experimental technique used to 
determine the form factor values. The KTeV detector, 
Monte Carlo simulation (MC), and data selection are 
described in Section IIVI Section explains the fitting 
procedure to extract the form factors and Section IVll dc- 
scribes the systematic uncertainties in the analysis. Fi- 
nally, in Section IVlII we present results for the form fac- 
tors and for the phase space integrals. 
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FIG. 1: Feynman diagram for the decay Kl — » 1X l^v. The 
hashed circle indicates the hadronic vertex that depends on 
decay form factors. 



II. FORM FACTOR PHENOMENOLOGY 



In the Standard Model, the Born-level matrix clement 
for the Kl — * n ± £ T i> decay modes (Fig. ^| is [32j 



M 



G 



-FV us [f + (t)(P K 

+f-(t)(Pi<- 



KP 7r )^7p( 1 + 75)^ 

PnYuil^ + lsjuAi 



(1) 



where Gf is the Fermi constant and /+ and /_ are the 
vector form factors. Time-reversal invariance guaranties 
that the form factors are real, while local creation of the 
lepton pair requires that they be functions only of the 
square of the four momentum transfer to the leptons: 



= (Pi + Puf = (Pk - P.) 2 = t 



The above equality holds for Born-level kinematics. For 
radiative Kl — > ir Pi/j events, we assume that the form 
factors depend only on the virtuality of the exchange bo- 
son at the hadronic vertex |33j. 

/+ and /_ are not the only possible pair of form factors 
that can describe semileptonic kaon decays. More recent 
experiments have chosen to use fa instead of /_ : 



fo(t) = /+(*) 



t 



Ml -Ml 



/-(*), 



(3) 



where Mk and M n are the kaon and pion mass. /+ 
and /o correspond to the 1~ (vector) and + (scalar) ex- 
change amplitudes, respectively. For the Born-level ma- 
trix element, /o is multiplied by the lepton mass Mf, for 
Kl — > ir e^i/, it is negligible compared to the contribu- 
tion of /+ . 

Note that for the /+(£) and fo(t) form factors, /+(0) = 
/o(0); otherwise f-(t) is divergent for t — ► 0. To simplify 
the notation, we use normalized form factors, /+o(£) = 
/+,o(*)//+(0). 

Historically, the normalized form factors f+(t) and 
fo(t) were usually analyzed assuming a linear t depen- 



U(t) = i + x. 

f (t) = 1 + A 



t 

~Ml 

Ml 



(4) 



This linear parametrization gave an adequate description 
of the experimental data. With a larger data sample and 
better control of systematic uncertainties, we study the 
second order term in /+ [34|: 



t 



1 



t z 



f +i t) = l + \>-, + ~\»-z 



Mi 



Mi 



(5) 



The presence of higher order terms is motivated by the 
pole model, in which the t dependence of /+(£) and fa(t) 
is described by exchange of the lightest vector and scalar 
K* mesons, respectively: 



hit) = 



M 2 
M 2 - 
M 2 



(6) 



Mi 



t 



Pole models expect M v rs 892 MeV, the mass of the 
lightest vector strange-meson, and M s > M v . The vec- 
tor meson dominance picture is supported by the recent 
results on r — ► Kttv decays showing a large enhancement 
of the K*(892) channel [aHOH 

From the Born-level matrix element (Eq. QJ, the 
Kl — > 7t ± ^ =f j/ decay rate is 



F 



M3 = ^fSEw(l + &) \Vu S \ 2 \fl(0)\I e K 



(7) 



Sew UJl is the short-distance radiative correction that is 
(2) the same for both modes and (l + S e K ) [l2,[lij are mode- 



dependent, long-distance radiative corrections. The di- 
mensionless decay phase space integrals, 1^, depend on 
the decay form factors |3| : 



I 



K 



dt j^g- A 3 / 2 F(tJ+(t)J (t)) , 



(8) 



where 

F(t,f + (t), f (t)) 



M 2 



M? 



f 3M 2 (M 2 K -M 2 Y f . 

I+[t)+ (2t + M 2 )X Io[t) 



(9) 



and 



\ = t z 



M 



K 



Ml - 2tM 2 K - 2tM 2 - 2M 2 K M 2 . (10) 



III. 



FORM FACTOR EXTRACTION 
PROCEDURE 



For a semileptonic decay detected by KTeV, the kaon 
direction of flight is well measured, but the kaon energy is 



not uniquely determined. Using the reconstructed miss- 
ing transverse momentum squared associated with the 
neutrino (p\ „) and the invariant mass of the pion-lepton 
system (Mrf), one can uniquely determine the square 
of the neutrino longitudinal momentum in the kaon rest 
frame: 



p* 2 



{Ml 



M^f 



&M K 



Pa 



(11) 



The sign ambiguity for P,f v leads to a twofold ambiguity 
for the parent kaon energy, Ek , which in turn gives rise 
to two different values of t. 

Most previous fixed target experiments encountered 
the same twofold kaon energy ambiguity. Typical solu- 
tions to this problem included selecting events with small 
|P|f J, picking the more probable solution, or using both 
solutions with weights proportional to the probability of 
each solution. In these approaches, the probability of a 
given solution depended on the form factors and the kaon 
energy spectrum, which was often not well understood; 
this coupling increased the systematic uncertainty. 

In this paper, we describe a technique to measure form 
factors that avoids the twofold Ek ambiguity. The key 
idea in this method is that instead of t, we use "trans- 
verse t," which is defined by substituting the measured 
transverse momenta for the full 3-momenta in Eq. [21 A 
similar approach was used in j 1 31 ] . 

We introduce two definitions of transverse t following 
the two definitions of t. We define "lepton transverse f\ 
ij_, by modifying (Pi + P v ) 2 : 



t^ = M$ + 2\p ±t „Up 2 ±4 + M%-2p± tV p± ! i 

= M} + 2|p ±jl/ |y^ + Mj +plj+pi, v - v\.« ■ 

(12) 
In a similar way, we define "pion transverse £," tj_, by 
modifying (P K - P^f: 



t\ = M 2 K + Ml - 2M K 0£ jW + Ml . (13) 

Note that from Eq. |3 the following inequality holds: 

*!<*<*!• (14) 

To simplify the notation, we use the symbol t± if the 
discussion refers to either t\ or t e ± ; the symbol tf is 
used to refer to both t\ and fj_. 

Since the Kl is a scalar particle, the decay kinemat- 
ics does not depend on the Kl direction. Therefore, the 
t± distribution can be unambiguously related to the t 
distribution and vice versa. The value of t± is invari- 
ant under kaon boosts, and thus the t±_ distribution does 
not depend on the kaon energy spectrum. The loss of 
longitudinal information does not dramatically increase 
the statistical uncertainty: a MC study shows that the 
£_i_-based measurement of the form factors leads to a 15% 
increase in uncertainty with respect to an ideal, i-based 
measurement. 
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FIG. 2: Plan view of the KTeV (E832) detector. The evac- 
uated decay volume ends with a thin vacuum window at 
Z = 159 m. The label "Csl" indicates the electromagnetic 
calorimeter. 



Each definition of tj_ can be used to measure the form 
factors. We will refer to the extraction of form factors 
using the t e ± or t\ variables as the t e ± - or ^-methods, 
respectively. The two methods have a different sensitiv- 
ity to systematic effects. For example, the ^-method is 
much less sensitive to radiative corrections, and therefore 
is used for the Kl — > 7r ± e =F ^ decay mode, while the t l ± - 
method leads to better statistical precision, and is used 
for the Kl —> it fjfiv decay mode. 

The form factor measurement starts with a determi- 
nation of the t± distribution in data. We use a Monte 
Carlo simulation to determine the detector acceptance 
as a function of t± and also to calculate radiative correc- 
tions. We extract the form factors by fitting the Monte 
Carlo simulation to data as will be described in SectionlVl 



IV. APPARATUS, SIMULATION, AND DATA 
SELECTION 

The KTeV experiment (Fig. |2J) and associated event 
reconstruction techniques have been described in detail 
elsewhere 15]. An 800 GeV/c proton beam striking a 
BeO target is used to produce two almost parallel neu- 
tral beams. A large vacuum decay region surrounded 
by photon veto detectors extends from 110 m to 159 m 
from the primary target. Following a thin vacuum win- 
dow at the end of the vacuum region is a drift chamber 
spectrometer, equipped with an analysis magnet which 
imparts a 0.41 GeV/c kick in the horizontal plane. Far- 
ther downstream, there is a trigger hodoscope, a pure Csl 
electromagnetic calorimeter, and a muon system consist- 
ing of two scintillator hodoscopes behind 4 m and 5 m of 
steel. 

A detailed Monte Carlo simulation is used to correct 
for the acceptance of the Kl — > 7r ± £ T i/ decay modes. The 
MC includes three main steps: event generation, propa- 
gation of particles through the detector, and simulation 
of the detector performance. A thorough discussion of 
this simulation program is given in pL51l , and more specif- 
ically for Kl — + Ti^F^v decays in [lfj. To include leading- 



order QED radiative corrections, we use the KLOR pro- 
gram [HI- 

For this analysis, we use a low intensity data sample 
collected in 1997, which is the same data sample used to 
measure T(K L -> tt ± ^v)/T{K l -> n ± e^i/) 10]. Con- 
sequently, most of the data selection requirements are 
similar to those used in the KTeV partial width ratio 
analysis. 

The reconstruction begins with the identification of 
two oppositely charged tracks coming from a single ver- 
tex, reconstructed between 123 and 158 m from the pri- 
mary target. To pass the event selection, one of the 
two tracks must be within 7 cm of a cluster in the Csl 
calorimeter; the second track is not required to have a 
cluster match. The event's missing transverse momen- 
tum squared associated with the neutrino is calculated 
with respect to the line connecting the primary target 
and the decay vertex. Both kaon energy solutions are 
required to fall between 40 and 120 GeV, where the kaon 
momentum spectrum is well measured. 

The decay Kl — ► 7r ± e =F ^ is identified using E/p, the 
energy reconstructed in the Csl calorimeter divided by 
the momentum measured in the spectrometer. One track 
is required to have E/p greater than 0.92 and the other 
track is required to have E/p less than 0.85. The first 
condition is satisfied for 99.8% of electrons and failed 
by 99.5% of pions. The second condition is satisfied 
for 99.1% of pions and failed by 99.93% of electrons. 
The resulting unambiguous identification of the pion and 



electron in the decay Kl 



TT^ 



v (swap probability 



< 0.01%) allows for the correct reconstruction of t^. 

The decay Kl — > it y^v is identified by requiring 
one of the tracks to have an energy deposit in the 
Csl less than 2 GeV J35J, and the other track to have 
E/p less than 0.85. This selection strongly suppresses 
Kl -^ tt e^v decays; however, since the Csl calorime- 
ter is only about 1/3 of a hadronic interaction length, 
the calorimeter selection is not very effective in reducing 
background from Kl — > 7r + 7r~7r° and Kl —> 7r + 7r~ . 

In both semileptonic decay modes, Kl — * tt + tt~tt° 

background is rejected using the variable fc_| o |l6j - 

For a genuine Kl — * 7r + 7r~7r° event, k^ o is propor- 
tional to the square of the ir° longitudinal momentum 
in the reference frame in which the sum of 7r + -7r~ mo- 
menta is orthogonal to the kaon direction. We require 
fc+_ < -0.006, rejecting 99.9% of K L -> vr+Tr^vr de- 
cays while retaining 99.5% of Kl — > w e^v and 97.9% of 
Kl —* vr n^v decays. Kl —* tt + tt~ background is sup- 
pressed by rejecting events in which the two-track in- 
variant mass, assuming both tracks are pions, is between 
0.488 and 0.505 GeV/c 2 . 

The main difference in the selection of the decay 
Kl — ► tt* yP'v compared to the partial width ratio analy- 
sis is driven by the need to distinguish between the pion 
and muon tracks in order to reconstruct t±. Since the 
calorimeter selection does not uniquely identify the pion 
and muon, we identify the muon by requiring the extrap- 
olation of a track to be near fired horizontal and ver- 
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FIG. 3: The invariant mass distribution (a) M^e for 
Kl — » n e^v events and (b) M^^ for Kl — > ir fj^v events. 
Data are shown as dots; MC is shown as a histogram. All 
analysis requirements have been applied. 



tical counters in the muon system. The quality of this 
"match" is measured by a xt variable that accounts for 
momentum-dependent multiple scattering in the steel in 
front of the muon system. We reject events with two 
matched tracks to suppress muon misidentification from 
tt — > [iv decays and from pion penetrations through the 
steel. The misidentification is suppressed further by re- 
quiring 0.5 m separation between the two tracks at the 
muon system (compared to the 0.15 m segmentation of 
the muon system counters). The efficiency of the match- 
ing requirement is increased by demanding that the muon 
track projection be at least 0.25 m inside the muon sys- 
tem outer aperture. To reduce the fraction of muons 
ranging out in the steel, we also require that the muon 
track momentum be greater than 10 GeV. These selec- 
tion criteria result in the unambiguous identification of 
the pion and muon in the decay Kl — ► tt fi^v, with a 
swap probability of less than 0.1%. 

There are two categories of background in this analy- 
sis: misidentified kaon decays and scattering background 
caused by a kaon scattering off a beamline element (ab- 
sorber or collimator). The misidentification background 
for the Kl — ► 7r ± e =F ^ decay mode is below 10~ 5 and is 
ignored. The background for the Kl — > 7r ± ^i T ^ decay 
mode is also small (10~ 4 ), but is not uniform in tf 
and must be subtracted. Kaon scattering off beamline 
elements, resulting in an incorrect measurement of p\, 
affects (0.11 ± 0.01)% of all events. Systematic uncer- 
tainties arising from the modeling of the background are 
discussed in Section IVI CI 

The main kinematic variables for the extraction of the 
semileptonic form factors are the transverse momenta- 
squared of all three particles, p\ v , p\ e and p± n , as well 
as the pion-lepton invariant mass, M T f. Distributions 
of M„£ for the Kl — > 7r ± e T ^ and Kl — ► 7r ± ^ =F ^ decay 
modes are shown in Figure |3] Distributions of p\ „, p\ i 
and p\ are shown in Figure 01 
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t±, which is used to construct the matrix A rg , for (a) tj_- 
method and (b) t^-method. The plot is based on the 
Kl — > ■k jX^v decay mode. 



V. FITTING 

We employ two different strategies to extract the form 
factors from the measured t± distribution. In the first 
method ("fit-A"), we use the MC to obtain the matrix 
A rg , which relates the generated f" = (Pk — P?r) 2 , g- 
index, to the reconstructed t^or t\, r-index. Note that 
detector resolution and radiative effects are included in 
this matrix. Figure [S] shows the distribution of V and 



(b) p ±tV for K L -> n fFv events, (c) p ±iB for K L -+tt e*u ^ that ig uged tQ calculate the matrix ^ |jjj 
events, (d) p±^ for Kl ™* tv fi^u events, (e) p 



for 

Kl — * 7r ± e =F i' events, and (f) p±,n fo r Kl — > tt* (jl^v events 
(d). Data are shown as dots, signal MC is shown as histogram. 
All analysis requirements have been applied. 



Using the matrix A rg , we construct the prediction for 
the binned distribution of the reconstructed t± : 



j^predA = Nq Y A 



F(t g ji red ( tg )jr d 



(**)) 



3 



i 9 



F(t g J^(t g )J MC (t g )) 



(15) 



Note that the M^i and p]_ distributions (Figs. 1314(1 
have different sensitivity to radiative effects and form 
factor values. For example, using the Kl — > 7r ± e =F i/ MC 
without radiative effects, the form factor obtained with 
the i^-method gives good data-MC agreement in p\ n , 
but poor agreement in the M xe and p\ e distributions. 
Conversely, the form factor obtained with the t^_-method 
gives good data-MC agreement in M ve andp^_ e , but poor 
agreement in p±„. The excellent data-MC agreement 
shown in all four distributions (M ve , p\ , p\ and p\ 
in Figs. 1314(1 is possible only if radiative effects are in- 
cluded in the MC. 



The measurement of the semileptonic form factors ben- 
efits greatly from the excellent resolution of the KTeV 
spectrometer (< 1%), and also from the extensive cali- 
bration performed for the measurement of e'/e [l5j . The 
resulting t± resolution is better than 1.5% on average. 



Here, F(t, /+(£), /o(£)) is the Born-level t distribution 
function (Eq. |§J) which depends on the form factors [33 ; 



/+ o are the form factors used in MC and /? r g are the 
form factors floated in the fit; t g are the bin centers and 
iVo is the overall normalization, which is also floated in 
the fit. For /^g = /+o, the prediction AP^* is simply 
the reconstructed t± distribution from MC. For different 
n""o , the prediction function (Eq. 115(1 is equivalent to 
"reweighting" the MC to the new values of the form fac- 
tors. The actual fit is the \ 2 minimization between the 
binned t± distribution in data and the prediction A/j^f 

In the second method ("fit-B"), we follow a common 
approach in which we fit the data t± distribution, cor- 
rected for acceptance and radiative effects, to the Born- 
level prediction. First, we define the acceptance At, in 
a given bin of t± b as the ratio of the number of recon- 
structed events to the number of generated events in the 
bin. Next, we define the radiative correction (1 + <5b) for 
each bin as the ratio of the number of generated events for 
the nominal MC (including radiative effects from KLOR) 



to the number of generated events for a separate Born- 
level MC; these MC samples are normalized using the 
global radiative correction (1 + 6^-). After correcting for 
acceptance and radiative effects, the number of corrected 
data events for a given £j_ bin is given by 



ArBorn 

iv ±,b 



N rec 



A b (l+6 b ) 



(16) 



where N^ is the number of reconstructed events. Both 
acceptance and radiative corrections depend on the form 
factors and are therefore determined in an iterative pro- 
cedure starting with the PDG values. Since this form 
factor dependence is relatively weak, only one iteration 
is required for convergence. 

To determine the prediction for the Born-level t± 
distribution, we start with a quadratic parametriza- 
tion of the form factors (Eq. |3J). We square the ma- 
trix element (Eq. QJ and collect terms proportional to 
X' + , X' + , A", A" , A ,A , A ',A ' , and cross terms. As a 
result, \A4\ 2 may be expressed as a second order function 
of A + , A , A + , A : 
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FIG. 6: Illustration of the influence of the form factors A+ = 
0.03 and Ao = 0.02 on the t ± and M w)1 distributions. The 
figure shows 0.03 x I~^ b /T±^ (solid lines marked as "A+") 

and 0.02 x T^ b /X±^ (dashed lines marked as "Ao") (a) as a 
function of t\ for decay Kl — * it e^u, (b) as a function of 
t^_ for decay Kl — ► it /X^u, and (c) as a function of M„ M for 
decay Kl — > tv fi^u. 



\M\ 2 =l+X' + 2+ + A Z° 



(17) 



where X are functions of the Dalitz plot variables. We 
generate K L — > it e^v and K L — > it fJ^v decays accord- 
ing to three-body phase space kinematics and numeri- 
cally integrate over I for events falling into the same t± 
bin. In this way, the prediction for the binned t± distri- 
bution can be expressed as 



Ki-predB 
JV ±,b 



(A + ,A , ...) — l± t b + X + l ± b + \ 2x h 



(18) 



where Z_L,t> are the results of the integration explained 
above. For the pole model fit, we first Taylor expand the 
pole parametrization (Eq. |SJ) up to second order in t to 
obtain 



A' - Ml and A"-2 M " 



(19) 



(and similar expressions for A and Aq as a function of 
Ms), and then use the prediction given by Eg. 1181 

The ratios X+T~J_ b /I±,b and XoXj_ b/-^-L,6 govern the in- 
fluence of the A+ and Ao parameters on the t^ distribu- 
tions. These ratios are illustrated in Figure for typical 
values of A+ and Ao- In a similar manner, we determine 
the influence of A+ and Ao on the M I(I distribution for the 
ir^pO^v decay mode (Figure El(c)). For the decay 
7T e^v : A+ = 0.03 corresponds to a linear increase 
of the decay rate as a function of t £ ± ; for tj_/M T 2 = 6, 
this increase is about 40%. For K L — > tt^/x^i;, the effect 
of Ao is non-linear and is more pronounced for large ij_ 
and small m 2 ,^. 

The advantage of fit-A is that there is no need for 
an iterative determination of the acceptance since the 
method is almost independent of the form factors used 
in the MC. The matrix A ro also utilizes the detector 



K L 

K L 



resolution information, which leads to better statisti- 
cal precision compared to fit-B (by about 10%). On 
the other hand, the advantage of fit-B is an easy gen- 
eralization to 2-dimcnsional fits, in particular to fit in 
(Mni,t±) space. Although fit-B has less precision than 
fit-A for 1-dimensional fits, the precision is improved in 
2-dimcnsional fits by using the different shapes of the Ao 
and A + distributions as a function of m 2 ^ (Fig. 0(c)). 

Both fitting methods are checked using MC samples 
generated with different values of form factors. The es- 
timation of the statistical uncertainty is verified by gen- 
erating 100 independent MC samples and studying the 
distribution of the fitted form factors. For the one di- 
mensional fits in t±, the form factors for K L —> jre^i/ 
and K L — > ir^p^v are extracted using both fitting tech- 
niques; these measurements are consistent within the un- 
corrected statistical uncertainty. 



VI. SYSTEMATIC UNCERTAINTIES 

The sources of systematic uncertainty for this analy- 
sis are the same as for the partial width ratio measure- 
ments 16], but their impact on the result is different. 
The form factor measurement does not depend on uni- 
form losses that do not bias kinematic distributions. For 
example, losses caused by hadronic interactions inside the 
spectrometer have virtually no influence. On the other 
hand, effects that cause misreconstruction of event kine- 
matics, such as calibration of the analysis magnet kick, 
are more significant. 

Table [I] summarizes the systematic uncertainties for 
the form factor measurements. For each decay mode and 
for each type of parametrization, the uncertainties are 



evaluated for the results obtained using both the t\ and 
ij_ methods. For the nominal KTeV result, we select the 
method which leads to the smallest total error and use the 
other method as a crosscheck. The nominal form factor 
result for Kl — * n^e^v is based on a 1-dimensional fit 
to the £j_ distribution; the nominal form factor result 
for Kl — > 7r ± /i =F zy is based on a 2-dimensional fit to the 
(-M^^ry distribution. 

The systematic uncertainty in the form factor resulting 
from a given effect is evaluated by turning this effect off 
in the MC simulation and determining the correspond- 
ing change in the form factors. The systematic uncer- 
tainty in the form factor is taken as this change times 
the relative uncertainty of the systematic effect. For ex- 
ample, turning off muon scattering in the steel changes 
Ao by 1.5 x 10~ 3 ; since we model the effects of scatter- 
ing with 20% precision fSec. IVIB 21) . the corresponding 
uncertainty on Ao is 0.3 x 10~ 3 . 

The following sections discuss systematic uncertainties 
in the order given in Table [I] 



A. Calibration &: Alignment 

1. Drift chamber alignment 

The dimensions of the drift chambers are known to 
20 /im from optical survey. The chamber locations along 
the kaon beam (^-direction) were monitored during the 
run and are known to 100 /im. In the transverse direc- 
tions, the chambers are aligned in-situ using dedicated 
muon runs with the analysis magnet turned off. The 
transverse position of the primary target is known to 
50 /xm precision by projecting the total momentum of 
Kl — ► 7r + 7r~ events to Z = 0. 

The alignment in the horizontal direction is most im- 
portant because it affects the track momentum measure- 
ment. The residual uncertainty in the alignment pro- 
cedure is 20 /im. The systematic uncertainty reported 
in Tabled is based on the sum in quadrature of the form 
factor changes if DC1 is shifted horizontally by 20 /xm, ro- 
tated in the XY plane by 20 /xrad, has non-orthogonality 
between vertical and horizontal wires of 50 /xrad, or is 
shifted in Z by 100 /xm. 



2. Magnet Kick Calibration 

The calibration of the analysis magnet transverse kick 
is performed with Kl —^ tt + it~ decays. The 30 keV un- 
certainty in the PDG value of Mk leads to a 0.01% un- 
certainty in the transverse kick. 



B. Detector Simulation 

1. Multiple Scattering and Resolution 

As explained in |lq , the amount of detector material 
and the influence of electromagnetic interactions (multi- 
ple scattering, bremsstrahlung and 5-ray production) on 
the detector acceptance are described to 10%. For the 
systematic uncertainty arising from modeling of the drift 
chamber position resolution, we assume a 5% uncertainty 
for the Gaussian part of the chamber response and a 10% 
uncertainty for the non-Gaussian tail induced by the dis- 
crete statistics of the ionization process. This estimate 
is based on a data-MC comparison of the distance be- 
tween the track and reconstructed hits. The chamber 
inefficiency caused by early accidental activity as well as 
hadronic interactions are negligible for this analysis. 



2. Muon System Simulation 

Several effects must be simulated to reproduce the effi- 
ciency of the muon system. These effects are muon scat- 
tering in the steel, modeling of the cracks between the 
muon system counters, and pion penetrations through 
the steel. 

Multiple scattering in the steel is modeled using a 
parametrization of a GEANT simulation [ItJ • It is checked 
with muons from Kl ~^ it n^v decays, in which the 
pion is unambiguously identifie d by requiring a hadronic 
shower in the Csl calorimeter [33] • Figure shows the 
inefficiency of the muon match as a function of the muon 
momentum. From the distribution of distances between 
extrapolated muon tracks and hit muon counters, we con- 
clude that scattering in the steel is described to better 
than 20% precision. The size of the cracks between the 
muon system counters (~ 1 mm) is determined to 50% 
precision using the muon runs. 

The probability that a pion penetrates through the 
steel and leaves a signal in the muon system is mea- 
sured in data using pions from fully reconstructed 
Kl — ► 7r + 7T~7r° decays. This effect is included in the 



MC, and is checked for Kl 



~-ji t v decays by study- 



ing the fraction of events in which both tracks match 
to hit muon system counters. For data and MC we ob- 
serve matches for both tracks in 2.17% and 2.25% of the 
cases, respectively. According to the MC, events with two 
matches result from pion decays (~ 85%), pion penetra- 
tions (~ 10%) and accidental hits (~ 5%). This study 
shows that the combined effect of pion penetrations and 
decays is modeled to better than 10%. 



3. Cut Variation 

We have modified various selection criteria to check the 
stability of the measured form factors. The general strat- 
egy is to relax or tighten one requirement at a time while 



TABLE I: Systematic, statistical and total uncertainties for the form factors. 
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FIG. 7: Inefficiency of the muon hodoscope matching re- 
quirement as a function of the muon momentum measured 
with Kl — ► n fi^u decays. Arrow indicates the momentum 
requirement. The inefficiency increases with momentum be- 
cause of the xfi definition. 



retaining all other cuts. These requirements include the 
track separation requirement at the drift chambers and 
other fiducial and kinematic cuts. 

Tests leading to significant changes are added to the 
systematic uncertainty. For example, an extra require- 
ment of no photon-like clusters reconstructed in the Csl 
calorimeter for the Kl ~^ it e*v decay mode leads to a 
(— 0.13±0.07) x 10 -3 change in A+, which is included as a 
systematic uncertainty. The systematic uncertainties in 
the form factors are estimated as the sum in quadrature 
of these variations. 



4- Beam Shape and Kaon Energy Spectrum 

The transverse shape of the kaon beam is defined by 
a system of collimating elements, the most important of 
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FIG. 8: Ring distributions for (a) Kl — > w e^v, and (b) 
Kl — > 71 iX^v decay mode. Dots show data, "MC sum" refers 
to signal plus background MC, and "MC Scat" shows the 
contribution from kaon scattering. The arrow indicates the 
analysis requirement. Events in the tail are mostly from kaon 
scattering. 



which is the defining collimator located 88 m from the 
primary target. While the geometry of the defining col- 
limator is well known, its exact alignment is harder to 
determine, leading to an uncertainty in the beam shape 
near the beam edge. The beam shape is studied using a 
"RING" variable, which is calculated using the X-Y ver- 
tex position projected along the kaon direction of flight to 
the calorimeter surface [l5|. Figure |S] compares data and 
MC RING distributions for the two semilcptonic modes. A 
RING value of cm 2 corresponds to the beam center and 
a RING value of ~ 100 cm 2 corresponds to the beam edge; 
the beam size at the calorimeter is about 10 x 10 cm 2 . 

In this analysis, we select events with RING < 121 cm 2 . 
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FIG. 9: Distributions for data, background MC and sum of 
background and signal MC for (a) t\/M„ in Kl — ' " 0=F ' 
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decay mode, (b) t±/M^ in Kl — > IK e^v decay mode, (c) 
t ± /M% in K L -> t&^v decay mode, and (d) i±/M% in 
A"l — > ir^^v decay mode. 
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FIG. 10: Radiative correction 1 + 5 = (Trad/cBorn, calculated 
as a function of (a) t e ± /M% and t ± /M% for Kl — > n^e^f, and 
(b) i x /M^ and tl/M* for Al -> it^yFv. 



Figure ED shows rj_ and i x distributions for the data, 
signal, and background MC. The systematic uncertainty 
in the form factors is assigned as 100% of the change 
if the misreconstruction background is not subtracted 
from the data. 



Larger RING values result mostly from kaon scattering in 
the defining collimator. To check if the beam edge affects 
the form factor measurement, we repeat the analysis with 
RING < 70 cm 2 and also with RING < 200 cm 2 . There is 
no statistically significant change in the form factors, so 
no systematic uncertainty is assigned. 

For kaon energies between 40 and 120 GeV, the kaon 
momentum spectrum is known to 0.01% per GeV [16J |. 
To estimate the impact of this uncertainty on the form 
factor measurement, we modify the MC energy spectrum 
linearly within this uncertainty. 



C. Background 

The fraction of events in which the parent kaon has 
scattered in the defining collimator, as well as the shape 
of the resulting p\ K distribution, are measured with 
fully reconstructed Kl — ► ir + ir~Tr° decays. The descrip- 
tion of scattering is checked using events reconstructed 
in the sidebands of the ring number distribution (Fig.|HJ); 
this comparison shows that the level of scattering is cor- 
rect to better than 10%. 

The misreconstruction background in the 
Kl — > tv fjj^v decay mode is from Kl — > 7r ± e^^, 
Kl -^ 7r + 7r~7r and Kl — ► tt + tt~ decays. This back- 
ground is determined using the MC normalized with 
our branching fractions [lq . and is then subtracted. 
While this background is small, it is not uniform in t±. 



D. Fitting 

We study the systematic uncertainty arising in the fit- 
ting procedure by comparing fit-A with fit-B and by vary- 
ing bin sizes. The results obtained with the fit-A and 
fit-B agree with each other within the uncorrelated sta- 
tistical uncertainty. 

To evaluate the systematic uncertainty caused by 
changes in binning, we start with the two dimensional 
fit-B and vary the number of bins. The number of M^i 
bins is varied between 1 and 130, and the number of tx 
bins is varied between 1 and 80. Altogether, we consider 
about 60 fits between these extremes. Using the devia- 
tion of each fit from the nominal fit, and the uncorrelated 
statistical uncertainties, we construct a x 2 irl , summing 
over all 60 fits. We evaluate this Xbin f° r a ^ n ^ parain- 
eterizations. For example, A+ determined in the linear 
model fit to Kl — ► 7r ± e =F ^ has y^ in /dof = 71.5/61, and 
Ao determined in the linear model fit to Kl 
has XuJdof = 85.9/60. 



7T jJL^l 



To account for the observed non-statistical change 
of the fit results with binning changes, indicated by 
Xbin/dof > 1, we assign an additional systematic error. 
The size of this uncertainty is determined such that ad- 
dition of this error (in quadrature) to the uncorrelated 
statistical uncertainty of each fit leads to Xbin/ do f = 1. 
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E. Radiative corrections 

Figure 1101 shows the influence of radiative correc- 
tions on the tf distributions for the Kl — > tr^H^v de- 
cay modes. Using the £ ± -method, radiative correc- 
tions are significantly larger for Kl — > n e^v than for 
K L — > 7r ± ^ T ^; using the i ± -method, radiative correc- 
tions are similar in both decay modes. 

For the Kl — > -k e^v decay mode, variables that are 
based on lepton kinematics have much larger radiative 
corrections than variables based on pion kinematics. 
Consequently, the i ± -method has much greater sensitiv- 
ity to radiative corrections than the t ± -method. For the 
t^-method, the difference in A+ measured with the nom- 
inal MC compared to a MC without radiative effects is 
8 x 10 -3 ; for the i ± -method, the corresponding difference 
is only 1 x 10~ 3 . Using the nominal MC, the agreement 
between the form factor values obtained using the v\ and 
t\ methods, AA+ < 0.5 x 10~ 3 at 68% c.l. fSec. lVllCt . 
shows that radiative corrections for the Kl — > 7r ± e =F ^ de- 
cay mode are described to better than 10%. 



For the K L 



~-{i T v decay mode, a sensitive check of 



radiative corrections is achieved by comparing Ao mea- 
sured from the t e ± and M I(I distributions in separate 
1-dimensional fits. For our MC without radiative ef- 
fects, the two fits return Ao values that are different 
by AA = (6.5 ± 1.3) x 10~ 3 , while the nominal MC 
gives good agreement: AA = (—0.1 ± 1.3) x 10~ 3 . 
This agreement shows that radiative corrections for the 
Kl — > 7r ± /i T ^ decay mode are described to better than 
20% precision. 
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FIG. 11: la correlation contours for the (a) linear fit between 
Ao and A+, (b) quadratic fit between Ao and A+, (c) quadratic 
fit between Ao and A" , and (d) quadratic fit between A+ and 
A' + . The dotted curve is for Kl — > tv e^v only, the dashed 
curve is for Kl — » 7r ^v only, and the solid curve is for the 
average of the two decay modes. The open and closed circle 
in (d) represents A' + , A" obtained from the Taylor expansion 
of the pole parametrization (Eg. 1191 around Mk* and around 
M v , respectively. 



F. Summary of Systematic Uncertainties 

The total systematic uncertainty in the form fac- 
tors is calculated as the sum of systematic uncertain- 
ties in quadrature (Table [IJ. For A+ ineasured in the 
Kl —* 7T e^v decay mode, there are comparable contri- 
butions to the uncertainty from magnet kick calibration, 
cut variation, and multiple scattering and resolution ef- 
fects. For Ao measured in the Kl — ► 7r /j^v decay mode, 
the largest systematic uncertainty is from modeling ra- 
diative corrections in the MC; this uncertainty is lim- 
ited by the statistical uncertainty of the test described in 
Sec. EH] 



VII. RESULTS 



Form Factor results 



After all selection cuts and background subtraction, we 
have 1950543 K L -> n ± e T iy and 1535951 K L -> i^^ 
decays, respectively. Using these data, we extract f+(t) 
in the Kl — > tt e^v decay mode, and both /+(£) and 
fo(t) in the Kl — * x fj^v decay mode. For the nom- 



TABLE II: Form factor results for both Kl — ► ^t^v decay 
modes and the average measured for Linear, Quadratic and 
Pole mode parametrization. The uncertainties correspond to 
the total uncertainties; the breakdown into statistical and 
systematic uncertainty is given in Table [I] The correlation 
coefficients between the form factors are given in Table Hill 
% 2 /ndf for individual modes correspond to statistical uncer- 
tainty only; % 2 / n df for the average includes uncorrelated sys- 
tematic uncertainty. See text for explanation of the averages. 





Kl — » 'K ± e^v\ Kl — > w^fl^v 


Average 


Linear Parametrization 
(xlO -3 ) 


A+ 
Ao 


28.32 ± 0.57 


27.45 ±1.08 
16.57 ±1.25 


28.13 ±0.51 
16.35 ±1.21 


X'fndf 


81.0/65 


240.4/236 


0.5/1 


Quadratic Parametrization 
(xl0~ 3 ) 


A' + 
A'l 
Ao 


21.67 ±1.99 

2.87 ±0.78 


17.03 ±3.65 
4.43 ± 1.49 
12.81 ±1.83 


20.64 ± 1.75 
3.20 ±0.69 
13.72 ±1.31 


X'fndf 


62.2/64 


230.7/235 


1.5/2 


Pole Parametrization 

(MeV/c 2 ) 


M v 

Ms 


881.03 ±7.11 


889.19 ±16.20 
1167.14 ±42.00 


882.32 ± 6.54 
1173.80 ±39.47 


X Z /ndf 


66.3/65 


234.7/236 


0.2/1 
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inal /+(£) measurement in Kl — > 7r ± e T i/, we use a 1- 
dimensional fit to the t\ distribution performed with fit- 
A. The choice of t\ over t e ± is based on the much smaller 
sensitivity to radiative corrections (see discussion in Sec- 
tion EHI and Fig.HSJ. For the nominal /+(*) and f (t) 
measurement in Kl — ► it ^v, we use a 2-dimensional fit 
to the (Mn-^ijJ distribution performed with fit-B. The 
2-dimensional fit uses the different shape of /+(£) and 
/o(i) in M nfJ , (Fig. EJ), improving the statistical precision 
and reducing the correlation between the form factors 
compared to a 1-dimensional fit to the t± distribution. 

The form factor results for the two decay modes and 
for the three different parameterizations are given in Ta- 
ble [n] the correlation coefficients are given in Table IIIII 
For all measurements, the statistical and systematic un- 
certainties are comparable. The A% 2 = 1 contours for 
the linear and second order fits based on Kl — > 7r ± £ =F i/ 
decays are shown in Figure^2 The quality of all fits is ac- 
ceptable. On the other hand, for the Kl — > ir^e^v decay 
mode, fits to /+(£) using the quadratic and pole param- 
eterizations improve the % 2 by 18 and 15 units, respec- 
tively, compared to the linear parametrization. Taking 
the systematic uncertainty into account, the significance 
of the non- linear term in the t dependence of / + is about 
4cr. 

The KTeV result for f+(i) is consistent with a pole 
model. The fit returns a pole mass M v = (882.32 ± 
6.54) MeV, in fair agreement with the lightest vector K* 
mass (891.66 ±0.26 MeV). 

Lepton universality requires that the /+(£) form fac- 
tors be equal for the Kl — ► it e^v and Kl — ► tt fj^f 
decay modes. Therefore, we average results for the two 
decay modes, taking into account statistical and system- 
atic correlations among the parameters, as well as cor- 
relation of the systematic uncertainties between the de- 
cay modes. The resulting averages are given in the last 
column of Table [HI All parameterizations of /+ (t) are 
consistent for the two decay modes. Note that although 
the Kl — ► tt e^v decay mode is insensitive to /o(t), the 
correlation between A+ and Ao (see Fig. 1111) results in 
an average Ao that is different from Ao deterimined only 
from Kl 



n (X^v. 



B. Determination of Phase Space Integrals 

Using the KTeV average values of the form factors, we 
calculate the decay phase space integrals (Eq. [SJ . The 
results obtained for the three different parameterizations 
are given in Table llVl 



The quadratic parametrization for f+(t) leads to in- 
tegrals that are about 1% lower than those for the lin- 
ear parametrization. For our extraction of \V US \ J4[, we 
use the quadratic parametrization since the second order 
term is observed with about 4cr significance. The resulting 



decay phase space integrals are 

I e K = 0.15350 ± 0.00044 ± 0.00095 
I& = 0.10165 ± 0.00039 ± 0.00070, 



(20) 



where the second error is an additional systematic uncer- 
tainty based on the difference between the quadratic and 
pole models. 



C. Crosschecks 

We have performed a series of crosschecks on the form 
factor measurements. A similar set of crosschecks is per- 
formed for the linear, quadratic and pole parameteriza- 
tion. Among these tests, the best statistical sensitivity 
is achieved for the linear fit to A + for Kl —> n e^v, and 
for the linear fit to Ao with fixed A+ for Kl — ► tt \£^v. 
For simplicity, we report crosschecks using these linear 
paramctrizations . 



1. Consistency among Methods 

As discussed earlier, we determine form factors us- 
ing both the t\- and ^-methods. For A + measured in 
the Kl — * n e^u decay mode, the difference between 
the nominal ^-method and the ^-method is AA + = 
(+0.4 ± 0.3 s t a t) x 10 -3 , where the error is the uncorre- 
lated statistical uncertainty estimated using an ensemble 
of MC samples. Recall from Sec. IVIEI that without ra- 
diative effects in the MC, the t\ and t e ± methods disagree 
by about 20cr. For Ao measured in the Kl — > ir y^v de- 
cay mode, the difference between the nominal ^-method 
and the ^[-method, AA = (+0.3 ± 0.5 stat ) x 10~ 3 , is 
also consistent with zero. Similar agreement is observed 
for the quadratic and pole parameterizations. 

Other crosschecks include one- and two-dimensional 
fits as well as different form factor fitting techniques (fit- 
A and fit-B) explained in Section These crosschecks 
do not show any systematic biases. 

We have checked the Kl — > 7r ± e =F i/ result using an in- 
dependent data sample collected with x 10 higher beam 
intensity. The difference between this analysis and the 
nominal one is AA+ = (—0.3 ± 0.5 s ^ a ^) x 10~ 3 . 

The muon identification in the Kl — * tt fj^iy decay 
mode is checked with an analysis in which we do not 
use the muon system to identify the muon, but instead 
identify the pion by requiring a hadronic shower in the 
Csl, which occurs for about 60% of events. The dif- 
ference between this analysis (which identifies the pion) 
and the nominal analysis (which identifies the muon) is 
AA = (-2.4+ 1.5) x 10~ 3 . Here, the error is the un- 
corrected statistical uncertainty (0.8 x 10~ 3 ), combined 
with the additional systematic uncertainty arising from 
the momentum dependent efficiency of the pion identifi- 
cation requirement (1.3 x 10 -3 ). 
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A^, 



ty e^v 



1T fji^v, correlation 
quadratic, and d) pole parametrizations. For the average of the two Kl — > tt t^v decay 



using quadratic parametrization. For Kl 



TABLE III: a) Total correlation coefficients for Kl 

coefficients are shown for b) linear, c 

modes, correlation coefficients are shown for e) linear, f) quadratic, and g) pole parametrizations 



a) 

at 



AV 



-0.97 



b) 



Ao 



-0.38 



c) 



A^ 

a; 



Ao A + 



0.65 
-0.75 



-0.96 



d) 



M v 



Ms 



-0.46 



e) 



Ao 



-0.36 



f) 



A'_ 



Ap A + 



0.34 
-0.44 



0.97 



g) 



M v 



Ms 



-0.40 



TABLE IV: Phase space integrals for Kl 
based on results in Tables ITT1 and ITTT1 

T e 



tv t^v decays, 



1 K 



Linear model 0.15507 ± 0.00027 0.10294 ± 0.00026 

Quadratic model 0.15350 ± 0.00044 0.10165 ± 0.00039 

Pole model 0.15445 ± 0.00023 0.10235 ± 0.00022 



2. Consistency among Data Subsets 

The stability of the results is tested by dividing the 
data into subsamples under a variety of criteria. These 
include decay vertex, minimum track momentum, track 
separation at the Csl calorimeter, pion-lepton mass, neu- 
trino direction, magnet polarity and lepton charge. All 
of these checks agree within the uncorrelated statistical 
uncertainty. A few particularly interesting checks are dis- 
cussed below. 

To verify that the tj^-methods are not biased by the 
ambiguity in the kaon energy, we divide the data into 
two nearly equal subsamples based on the direction of 
the neutrino in the kaon rest frame: 



I cos 6* 



\ 



p* 2 
^11." 



(P* 2 



II." 



PlvY 



(21) 



For cos(8*) = 0, the neutrino is emitted perpendicular 
to the kaon flight direction, and there is only one kaon 
energy solution; for |cos(6**)| = 1, the two kaon energy 
solutions have the maximum difference. The two sub- 
samples for this test are selected by | cos(0*)| < 0.4, for 
which the two kaon energy solutions differ on average by 
18%, and by | cos(0*)| > 0.4, for which the two kaon en- 
ergy solutions differ by 40%. For these two K L — > ir ± e T v 
samples, we measure A + = (28.0 ± 0.5 s ^ a t) x 10~ 3 
and A + = (28.9 ± 0.7 s t a t) x 10~ 3 , respectively, show- 
ing good agreement. For K L ~ * n^fj^v, there is also 
good agreement: A o (|cos0*| < 0.4)— A (| cos6>*| > 0.4) = 
(l-4±l-8 stat )xl0- 3 . 

We check the horizontal alignment of the spectrom- 
eter by dividing the K L —* it e^v data into four sub- 
samples based on the lepton charge and the analyz- 
ing magnet polarity. This procedure separates the data 
into classes with tracks bending in different directions 
in the horizontal plane. We find that for the nominal 
alignment, the four measurements of A+ are consistent 
( x 2 /ndf = 2.95/3); if the first drift chamber (DC1) is 
shifted by 50 /im in x-direction, the agreement is signif- 
icantly worse (x 2 /ndf — 30.90/3). Good agreement is 



also observed for the same test performed for Ao mea- 
sured in the Kl —> 7r ± /i =F ^ decay mode: \ 2 l n df — 1.8/3. 



3. Lepton Universality 

Lepton universality implies that the coupling constant 
Gf and the short-distance radiative correction Sew are 



the same for Kl — > 7r ± e =F ^ and Kl 



tx fj,^l 



Taking 



the ratio of Eq. 0for these two modes, we obtain a pre- 
diction for the partial width ratio: 



1 



[^K fJ ,3/^Ke3]p TC d - Y 



■% 



'K 



JK 

1 K 



(22) 



where the ratio of radiative corrections is calculated to be 
( 1 +Sk)/^ + S c k ) = 1.0058±0.0010|I3. The ratios of the 
phase space integrals calculated for the linear, quadratic 
and pole models are: 



1 K 
T e 



0.6639 ± 0.0017 Linear parametrization 

0.6622 ± 0.0017 Quadratic parametrization 



0.6627 ±0.0015 



Following the prescription 
quadratic parametrization, 
error based on the pole 

J&Ak = °- 6622 ± °- 0018 



Pole parametrization . 

(23) 

of Sec IVlTBl we use the 

and include and additional 

parametrization to obtain 

Using this ratio of the in- 



tegrals, and the KTeV value of r i ^ 3 /r is:e 3=0.6640 ± 
0.0014 ± 0.0022 [3, we obtain 



^Kfj.3 /Tk 



r.'i 



[^K^/^Ke3]p Te ^ 



0.9969 ±0.0048 



(24) 



which shows that the KTeV form factors are consistent 
with the KTeV partial width ratio measurement. 



D. Comparison with other Form Factor 
Measurements 



Previous experiments analyzed their data using a lin- 
ear approximation for /+(£) and fo(t)- To compare 
with these results, we also use the linear parametriza- 
tion. The comparison of results for A + measured in the 
Kl —> "K €*v decay mode is shown in Fig. El the A+ 
comparison in Kl — > 7r ± /i =F j/ is shown in Fig. 1131 and 
the Ao comparison in Kl — > ir^yAv is shown in Fig. 1141 
To simplify the comparison of Ao, all measurements are 
adjusted to the same value of A + = 0.030 using the re- 
ported dependences |39| . 
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Birulev 81 
CPLEAR 00 h 
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FIG. 12: A+ measured in Kl — > 7r e^v decay mode by KTeV 
and by experiments used in the PDG evaluation (Chien 71 
Hi|, Bisi 7 1 JTsl . Buchanan 75 HJ, Blumenthal 75 El. 
Gjesdal 76 [2l~Hill 78 O, Engler 78 0, Birulev 81 H|, 
CPLEAR 00~H)- The PDG fit (excluding KTeV) is shown 
by the vertical lines. 



Donaldson 74 

Buchanan 75 

i 

Hill 79 i 



Cho80 



Birulev 81 



KTEV 



0.02 0.03 0.04 



0.05 



FIG. 13: A+ measured in Kl — > vr \H^v decay mode by KTeV 
and by experiments used in the PDG evaluation (Donaldson 
74pf , Buchanan 75 [H, Hill 79 |27|, Cho 80 HI, Birulev 81 
|25|1. The PDG fit (excluding KTeV) is shown by the vertical 
lines. 



Sandweiss 73 

o 1 

Clark 77 
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X (for ?i + =0.030) 



FIG. 14: Ao measured by KTeV and by experiments used in 
the PDG evaluation (Sandweiss 73 [2<j, Clark 77 [H, Don- 
aldson 74 E|, Buchanan 75 |2JJ, Hill 79 ^, Cho 80 H|, 
Birulev 81 (Hi). PDG BR and KTeV BR are determined us- 
ing r/f M 3/Fif e 3 and lepton universality. KTeV Dalitz refers 
to the technique described in this paper. For comparison, all 
central values and uncertainties are adjusted to A+ = 0.030, 
using reported dependencies. The PDG fit, based on /i- 
polarization and Dalitz plot measurements (excluding KTeV 
result) is shown by the vertical lines. The curve represents 
a PDG-style ideogram constructed for all but the two KTeV 
results. 



The KTeV A+ result for K L -^-n^e^v is three 
times more precise than the PDG evaluation; for 
K L — > 7r ± /x =F ^, our result is almost five times more pre- 
cise. The new KTeV measurements of A + for both 
semileptonic decay modes are in agreement with most 
previous experiments and with the PDG average. The 
only significant disagreement is between the measure- 
ments of A-|- for K L — > 7T n^v performed by KTeV and 
Birulev 81 |2£|; the difference between these two mea- 
surements is more than 3<r (Fig 113ft . 

For Ao, the experimental situation is less clear 
(Fig. I14fl . The KTeV result, which is nearly five times 
more precise than the PDG evaluation, agrees with only 
two of the eight previous measurements. Figure 1141 
also shows a PDG-style ideogram calculated without the 
KTeV result. This ideogram indicates clustering of the 
measurements into two groups, one around Ao ~ 0.020 
( "middle Ao group") and another around Ao ~ 0.040 
("high A group"). 

Three different methods have been used to extract Ao: 
muon polarization, Dalitz plot density, and Tk^/^ Ke3 
partial width ratio together with the lepton universal- 
ity constraint (Eq. 1221 . The largest inconsistency is a 
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5<7 discrepancy between the two measurements based on 
muon polarization (Sandweiss 73 J2JJ and Clark 77 |30|'). 
Sandweiss 73 is the only result with a negative Ao- The 
KTeV result does not agree with either of these measure- 
ments. 

Among the Ao measurements based on the Dalitz plot 
analysis, KTcV agrees well with the highest precision 
measurement of Donaldson 74 26], but disagrees with 
the more recent measurements (Hill 79 J53, Cho 80 [28J . 
Birulev 81 25]). The largest inconsistency is a 4er dis- 
crepancy between KTeV and Cho 80 [23] ■ 

The evaluation of Ao from the KTeV measurement of 
^Kfiz/^KeZ agrees with the middle Ao group, while the 
same evaluation based on the PDG fit lands in the high 
Ao group. The PDG value for TK^/^KeS is to a large ex- 
tent driven by the measurement reported by Cho 80 |28l ] , 
which disagrees with the analogous KTeV measurement 
by more than 3a. 

It is worth mentioning that the physics interpretations 
of the middle and high Ao groups are quite different. For 
the former, the form factors are consistent with a pole 
model (Aq < A+), and also agree with those measured 
for K J40j, suggesting that isospin symmetry breaking 
effects are small. For the high Ao group, where Ao > A + 
(2a), the simple pole model is disfavoured; the high Ao 
group also shows more than 3a significance of isospin 
symmetry breaking for fo(t). 



PDG average. KTeV measurements of /+(£) are in good 
agreement with those used in the PDG evaluation. 

The KTeV result for Ao, which is nearly five times more 
precise than the PDG average, disagrees with six of the 
eight measurements used in the PDG average. Despite 
the poor agreement between KTeV and previous results, 
the most precise measurement |2g used in the PDG eval- 
uation agrees with the KTeV measurement. 

The new KTeV result indicates the presence of a non- 
linear term in the t dependence of /+ with about Aa 
significance. This non-linear term is consistent with pole 
model expectations. The KTeV result for Ao is also con- 
sistent with the pole model prediction, and is in agree- 
ment with the value measured in K decays, suggesting 
that isospin breaking effects are small. 

Using the KTeV measurements of the semileptonic 
form factors, we have calculated the values of the phase 



space integrals for the Kl 



z e T v and Kl 



7T fj/^l 



decay modes. We find that inclusion of the non-linear 
term in the /+(£) expansion reduces the value of these 
integrals by about 1%. Lepton universality holds for the 
KTeV data with 0.5% precision, showing consistency be- 
tween the KTeV form factor measurements presented in 
this paper^ KTeV semileptonic partial width ratios re- 
ported in [lj| , and radiative corrections calculated in [13j . 



VIII. CONCLUSIONS 



In this paper, we have presented new measurements of 
the form factors in semileptonic Kl decays performed for 
Kl — » 7r e^v and Kl — ► ^ \i^v. The measurement of 
/+(£) is consistent for the two decay modes, and provides 
about a threefold increase in precision compared to the 
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